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Cross-sectional  scanning  tunneling  microscopy  of  Mn-doped  GaAs:  Theory  and  experiment 
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We  report  first-principles  calculations  of  the  energetics  and  simulated  scanning  tunneling  microscopy  (STM) 
images  for  Mn  dopants  near  the  GaAs  (110)  surface,  and  compare  the  results  with  cross-sectional  STM  images. 

The  Mn  configurations  considered  here  include  substitutionals,  interstitials,  and  complexes  of  substitutionals 
and  interstitials  in  the  first  three  layers  near  the  surface.  Based  on  detailed  comparisons  of  the  simulated  and 
experimental  images,  we  identify  three  types  of  Mn  configurations  imaged  at  the  surface:  (1)  single  Mn 
substitutionals,  (2)  pairs  of  Mn  substitutionals,  and  (3)  complexes  of  Mn  substitutionals  and  interstitials. 
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I.  INTRODUCTION 

GaAs  can  be  doped  with  Mn  to  form  a  dilute  magnetic 
semiconductor  with  a  Curie  temperature  as  high  as  160  K.1  3 
It  is  generally  accepted  that  holes  created  by  the  substitution 
of  Mn  for  Ga  mediate  the  ferromagnetic  interaction  between 
Mn  dopants  in  this  material.4'5  Naively,  each  substitutional 
Mn  is  expected  to  produce  one  hole;  therefore  the  nominal 
hole  concentration,  p,  should  be  equal  to  the  Mn  concentra¬ 
tion.  For  typical  Mn  fractions  of  5%,  this  implies  a  hole 
concentration  p  =  1 . 1  X  1 02 1  cm-3.  Measured  hole  concen¬ 
trations  of  as-grown  material  are  much  less  than  this,  typi¬ 
cally  by  a  factor  of  ~3.  The  source  of  this  compensation 
remains  somewhat  controversial,  having  been  variously  at¬ 
tributed  to  either  excess  As  in  the  form  of  antisites  and 
interstitials6-9  or  Mn  interstitials.10-13  Recent  experiments 
show  that  careful  annealing  near  the  growth  temperature  can 
significantly  enhance  the  conductivity  and  hole 
concentration.3 11  This,  as  well  as  recent  ion  channelling 
experiments,11  suggests  that  interstitial  Mn  is  the  more  likely 
source.  A  complete  picture  is  lacking,  however;  for  example, 
the  distribution  of  Mn  interstitials  in  as-grown  material  is 
still  unclear. 

Cross-sectional  scanning  tunneling  microscopy  (XSTM) 
is  an  effective  tool  for  addressing  this  issue  because  it  can 
image,  with  atomic  resolution,  the  structural  and  electronic 
configuration  of  impurities  and  defects  as  they  are  present  in 
the  bulk.  There  are  several  recent  XSTM  studies  of  Mn- 
doped  GaAs;13-16  these  studies  used  the  correlation  between 
transport  data,  level  of  Mn  doping,  and  observed  defect  den¬ 
sity  to  infer  the  location  (substitutional  versus  interstitial)  of 
Mn  dopants  and  As  antisites.  The  relative  abundance  of  As 
antisites  and  interstitial  Mn  deduced  by  these  studies  lacks  a 
clear  consensus;  one  study  finds  that  antisites  are  not  present 
at  all,15  while  another  finds  that  both  antisites  and  interstitial 
Mn  are  present,  and  compensate  substitutional  Mn.13 

To  investigate  the  nature  of  dopants  and  defects  in  this 
material,  we  performed  high-resolution  XSTM  measure¬ 
ments  on  a  (110)  cleavage  plane  of  Mn-doped  GaAs  and 
used  first-principles  calculations  to  interpret  the  images.  Spe¬ 
cifically,  we  used  density-functional  theory  to  simulate  STM 
images  for  a  number  of  Mn  configurations  near  the  GaAs 
(110)  surface.  Configurations  were  chosen  based  on  their  cal¬ 
culated  energetic  stability,  and  included  several  metastable 


configurations  that  would  probably  be  kinetically  stabilized. 
A  detailed  comparison  of  the  resulting  simulated  and  experi¬ 
mental  images  leads  us  to  the  following  conclusions:  (1)  iso¬ 
lated  substitutional  Mn  as  well  as  pairs  of  substitutional  Mn 
occur  with  roughly  comparable  frequency  in  bulk  Mn-doped 
GaAs  and  (2)  interstitial  Mn  is  typically  bound  in  complexes 
with  substitutional  Mn  in  several  different  configurations. 

II.  EXPERIMENTAL  XSTM  IMAGES 

The  Mn-doped  GaAs  sample  used  in  this  study  was 
grown  by  molecular-beam  epitaxy  (MBE)  using  well  estab¬ 
lished  methods  and  conditions.17  A  buffer  layer  of  high  qual¬ 
ity  undoped  GaAs  was  first  grown  on  an  n-type  GaAs  sub¬ 
strate.  This  was  followed  by  growth  of  a  100-nm-thick  GaAs 
buffer  layer  grown  at  250  °C  and  a  260-nm-thick  layer  of 
Mn-doped  GaAs  at  the  same  temperature.  Growth  quality 
was  monitored  by  reflection  high-energy  electron  diffraction. 
The  Mn  concentration  was  determined  with  x-ray  diffraction 
to  be  ~  1%. 

This  sample  was  then  characterized  by  XSTM  measure¬ 
ments  in  a  multichamber  ultrahigh-vacuum  facility.  In  order 
to  obtain  an  atomically  abrupt  cleave  across  the  epilayer,  the 
samples  were  thinned  ex  situ  to  ss  200/xm.  After  being 
loaded  into  the  XSTM  chamber  (base  pressure  =£10-10 
Torr),  the  samples  were  scribed  and  cleaved  in  situ  to  expose 
a  (110)  surface  (perpendicular  to  the  MBE  growth  direction), 
as  previously  described.13  The  constant-current  (40  pA)  im¬ 
ages  shown  here  are  of  filled  electronic  states  measured  with 
a  sample  bias  of  —2.5  V. 

Figure  1  shows  the  results  of  filled-state  XSTM  measure¬ 
ments  on  this  sample.  The  scan  area  of  this  image  is  170 
X  170  A2  and  displays  ~  1500  surface  atoms.  Filled-state 
STM  images  of  GaAs  (110)  reveal  the  As  surface  sublattice; 
the  Ga  sublattice  evident  in  empty-state  images,  on  the  other 
hand,  is  not  directly  revealed  in  this  image.  Since  only  the 
top  As  atoms  contribute,  STM  of  the  GaAs  (110)  surface 
reveals  only  every  other  layer  in  the  (001)  direction. 

Figure  1  reveals  many  surface  features  not  normally  seen 
in  STM  of  the  GaAs  (110)  surface.  Several  occur  with 
enough  regularity  to  warrant  close  attention.  Feature  A  is  an 
elliptical  region  of  intensity  in  an  As  (001)  plane  with  major 
and  minor  axes  the  size  of  two  and  one  As  surface  atoms, 
respectively.  Feature  B  is  circular,  and  only  slightly  larger 
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FIG.  1.  Constant-current  filled-state  XSTM  image  of  a  (110) 
cleavage  plane  of  Mn-doped  GaAs  (001).  Five  reproducible  fea¬ 
tures  (A-E)  are  marked  on  the  image  and  drawn  schematically 
below. 

and  more  intense  than  the  contributions  from  the  surrounding 
As.  Feature  C  is  more  diffuse  and  broad  than  both  A  and  B, 
with  intensity  spread  over  two  pairs  of  As  atoms  in  neigh¬ 
boring  (001)  planes.  As  discussed  in  more  detail  below,  this 
is  what  one  would  expect  if  there  were  two  A  features  in 
neighboring  As  (001)  planes.  Feature  D  is  of  similar  shape  to 
that  of  C  but  is  more  asymmetric,  with  one  side  of  the  defect 
the  size  of  a  single  surface  As  atom.  There  does  not  appear  to 
be  a  preferred  crystalline  orientation  of  feature  D,  as  differ¬ 
ent  orientations  are  visible  in  Fig.  1.  Feature  E  is  a  dark 
region  close  to  a  surface  As  atom  displaced  in  the  (00T) 
direction;  there  is  also  an  overall  apparent  depression  in  the 
vicinity  of  this  feature,  characteristic  of  band  bending  that 
occurs  near  positively  charged  defects  on  a  p-typc  surface.19 

On  the  basis  of  the  experimental  data  alone  it  is  extremely 
difficult,  if  not  impossible,  to  determine  what  dopants  or  de¬ 
fects  produce  these  features.  As  we  show  in  this  work,  how¬ 
ever,  when  interpreted  with  first-principles  theory,  these  mea¬ 
surements  resolve  both  the  type  and  geometric  configuration 
of  Mn  dopants  in  GaAs. 

III.  THEORETICAL  APPROACH 
A.  Physical  and  electronic  structure 

We  modeled  the  GaAs  (110)  surface  in  a  supercell  slab 
geometry  consisting  of  five  layers  of  GaAs  in  a  4X4  surface 
unit  cell.  A  vacuum  region  of  13.4  A  was  used  for  all  the 
calculations.  In  such  a  cell  the  separation  between  a  dopant 
and  its  periodic  image  is  —16  A,  sufficiently  large  to  ensure 
that  dopants  interact  only  negligibly  with  their  periodic  im¬ 
ages.  The  bottom  layer  of  the  slab  was  passivated;  the  atomic 
positions  were  fully  relaxed  in  all  but  the  bottom  two  layers. 

The  wave  functions  and  charge  density  were  expanded  in 
a  plane-wave  basis  and  evaluated  using  an  ultrasoft  pseudo- 
potential  approach  as  implemented  in  the  VASP  code." 


Electron  correlations  were  treated  at  the  level  of  the  local 
spin-density  approximation  with  total  energies  converged  to 
1(T4  eV. 

A  single  wave  vector  was  used  to  sample  the  Brillouin 
zone  of  the  slab  supercells  and  a  plane-wave  cutoff  of  227.24 
eV  was  used  for  all  the  calculations.  For  supercells  contain¬ 
ing  Mn,  atoms  within  —7.0  A  of  the  Mn  site(s)  were  fully 
relaxed  until  the  total  energy  between  two  structural  configu¬ 
rations  changed  by  less  than  10~3  eV. 

B.  Simulated  STM  images 

Theoretical  STM  images  were  simulated  using  the 
method  of  Tersoff  and  Hainan n.25  In  this  approach,  the  cen¬ 
tral  quantity  is  the  local  density  of  states  (LDOS)  near  the 
surface.  The  LDOS  is  integrated  over  an  energy  range  deter¬ 
mined  by  the  experimental  bias  voltage;  contours  of  constant 
energy-integrated  LDOS  simulate  a  constant-current  image. 
For  filled  state  imaging,  the  upper  bound  of  this  energy  range 
is  the  Fermi  level.  The  Mn-doped  GaAs  sample  under  con¬ 
sideration  in  this  work  is  p  type  with  an  estimated  hole  den¬ 
sity  of  1019  cm~3,  so  that  the  Fermi  level  should  be  near  the 
valence-band  maximum  (VBM).  Since  we  do  not  know  the 
distribution  of  donors  and  acceptors  present  in  the  sample, 
this  information  alone  is  not  sufficient  to  determine  the  exact 
location  of  the  Fermi  level  for  our  sample.  However,  we 
have  found  that,  as  a  consequence  of  the  large  bias,  the  simu¬ 
lated  images  are  rather  insensitive  to  the  value  of  the  Fermi 
level  used  in  the  LDOS  integration,  and  so  we  choose  arbi¬ 
trarily  a  point  0.4  eV  above  the  VBM  as  the  upper  bound  for 
all  the  STM  simulations.  For  this  choice,  the  Fermi  level  is 
above  the  acceptor  level  of  substitutional  Mn  and  below  the 
donor  levels  of  interstitial  Mn.  Simulations  performed  with 
the  position  of  the  Fermi  level  shifted  by  ±0.3  eV  give  simi¬ 
lar  results.  The  location  of  the  VBM  for  each  type  of  dopant 
was  determined  by  inspection  of  the  total  DOS  relative  to 
that  projected  onto  the  dopant  site. 

As  a  measure  of  the  agreement  between  simulated  and 
measured  STM  images,  we  consider  mainly  the  overall  shape 
and  spatial  extent  of  features  of  interest.  The  defect-free 
(110)  surface  has  a  (1 10)  mirror-plane  symmetry,  so  it  is 
also  useful  to  note  whether  a  particular  dopant  configuration 
preserves,  even  approximately,  this  symmetry.  Finally,  we 
use  qualitative  arguments  to  judge  whether  the  number  of 
observed  features  are  consistent  with  our  total-energy  calcu¬ 
lations. 

IV.  SUBSTITUTIONAL  Mn 
A.  Simple  substitutional  dopants 

Substitutional  Mn  dopants  are  of  primary  importance  for 
the  magnetic  and  transport  properties  of  Mn-doped  GaAs, 
and  hence  we  discuss  the  simulated  STM  images  for  these 
types  of  dopants  first.  Since  STM  is  a  near-surface  sensitive 
probe,  we  have  only  considered  substitutional  Mn  dopants  in 
the  top  three  layers  near  the  (110)  surface,  as  shown  in  Fig. 
2(a). 

We  will  use  the  notation  s(n)  to  denote  a  substitutional 
Mn  in  layer  n  (with  the  topmost  surface  layer  defined  as  n 
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FIG.  2.  (Color  online)  Top  view  (upper)  and  side  view  (lower) 
showing  the  Mn  dopant  configurations  considered  in  this  work.  Ga, 
As,  and  Mn  atoms  are  shown  as  cyan,  yellow,  and  red  circles, 
respectively.  The  dashed  lines  separate  regions  containing  different 
types  of  dopants.  From  left  to  right:  (a)  substitutional  Mn  in  the  first 
three  layers,  (b)  interstitial  Mn  in  the  first  two  layers,  (c)  complexes 
of  substitutional  and  interstitial  Mn  in  layer  2,  (d)  complexes  of 
substitutional  Mn  in  layer  1  and  interstitial  Mn  in  layer  2,  (e)  com¬ 
plexes  of  substitutional  and  interstitial  Mn  both  in  layer  1.  For  the 
complexes  in  regions  (c),  (d),  and  (e),  several  alternative  locations 
of  the  interstitial  Mn  involved  in  the  complex  are  shown  in  a  lighter 
shade  of  red.  See  text  for  explanation  of  the  notation. 


=  1).  Figures  3(a),  3(b),  and  3(c)  show  the  simulated  images 
from  s(l),  s(2),  and  s(3),  respectively.  The  s(l)  and  s( 3) 
images  are  strikingly  similar,  for  the  following  reasons.  ,v ( 1 ) 
forms  bonds  with  two  As  atoms  that  are  nominally  in  the 
same  plane;  likewise  for  the  s( 3)  Mn  atom.  Thus  the  overall 
shape  of  the  dopant  image  consists  of  three  overlapping 
maxima,  centered  on  the  two  As  atoms  and  on  the  Mn  atom, 
giving  rise  to  an  elongated  feature  with  a  (IlO)  mirror-plane 
symmetry.  Hence,  both  of  these  sites  are  consistent  with  fea¬ 
ture  A  in  Fig.  1. 

.v ( 2 )  on  the  other  hand,  has  a  very  different  image.  ,v ( 2 ) 
forms  a  bond  with  an  As  atom  in  the  top  layer,  with  the  bond 
pointing  toward  the  surface.  The  perturbation  due  to  s(2)  is 
centered  on  the  surface  As  site  involved  in  the  bond.  This 
overlaps  with  the  contribution  from  the  Mn  dopant  itself,  and 
so  the  shape  is  more  circular  compared  to  ,v  ( 1 )  and  ,y  ( 3 ) . 
Thus  it  enhances  a  single  As  site  and  has  a  ( 1 10)  mirror- 
plane  symmetry.  These  agree  qualitatively  with  the  charac¬ 
teristics  of  feature  B  in  Fig.  1,  leading  us  to  identify  those 
features  as  arising  from  ,v ( 2 )  dopants. 


B.  Substitutional  Mn  pairs 

Since  there  will  always  be  some  fraction  of  neighboring 
Mn  that  occur  simply  by  chance,  we  consider  pairs  of  sub¬ 
stitutional  Mn  as  well.  Figures  3(d),  3(e),  and  3(f)  show 


FIG.  3.  Panels  (a),  (b),  and  (c)  show  simulated  STM  images  of 
isolated  Mn  substitutionals.  Panels  (d)  through  (h)  show  simulated 
images  for  Mn  substitutional  pairs.  Mn  positions  are  marked  with 
+ .  The  crystalline  orientation  is  the  same  as  in  Fig.  2. 

simulated  images  of  pairs  of  Mn  substitutionals  s(l)s(l), 
s(3)s(3),  and  s(3)s(l),  respectively.  For  the  s(3)s(l)  pair, 
we  note  that  even  though  the  two  Mn  involved  are  ~6.8  A 
apart,  and  in  different  atomic  planes,  because  of  the  pro¬ 
jected  view  of  the  STM  scan  the  appearance  of  this  pair  is 
quite  similar  to  those  of  the  ,y(  1  )s(  1 )  and  ,y ( 3 ) ,y ( 3 )  pairs.  All 
three  of  these  simulated  images  have  a  shape  comparable  to 
that  of  feature  C  in  Fig.  1.  Common  aspects  include  a  region 
of  low  intensity  between  the  two  dopants,  a  (TlO)  mirror- 
plane  symmetry  (the  experimental  data  only  approximately 
have  this  symmetry)  and  a  similar  spatial  extent  involving 
four  surface  As  atoms.  It  is  also  interesting  to  note  that  these 
pairs  have  simulated  images  that  are  essentially  a  linear  su¬ 
perposition  of  the  images  of  the  isolated  substitutionals 
shown  in  Figs.  3(a)  and  3(c).  This  suggests  that,  at  the  level 
of  detail  that  STM  reveals,  two  substitutional  Mn  do  not 
interact  very  strongly  when  they  are  next-nearest  neighbors. 

Figures  3(g)  and  3(h)  show  the  simulated  images  for 
nearest-neighbor  substitutional  pairs  ,y  ( 1 )  .y  ( 2 )  and  ,y  ( 3 )  s  ( 2 ) , 
respectively.  Overall,  the  simulated  images  of  these  pairs  are 
similar  to  what  one  would  obtain  from  a  linear  superposition 
of  the  images  of  the  isolated  defects.  However,  the  simulated 
image  for  .?( 1  )s(2)  clearly  reveals  that  the  contribution  from 
the  s(l)  dopant  is  now  asymmetric,  lacking  a  (110)  mirror 
plane.  This  indicates  that  substitutional  pairs  interact  more 
strongly  when  they  are  in  nearest-neighbor  configuration 
than  in  the  next-nearest-neighbor  configuration  described  in 
the  paragraph  above.  The  shape  of  these  simulated  images  is 
generally  consistent  with  that  of  feature  D  in  Fig.  1,  with  a 
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characteristic  “L”  shape  involving  three  As  atoms,  two  of 
which  are  in  the  same  As  (001)  plane.  We  have  not  found  any 
other  dopants  that  yield  such  a  favorable  comparison  to  the 
experimental  data,  leading  us  to  identify  feature  D  as 
,v ( 1 ) ,v ( 2 )  and  ,v ( 3 ) ,v ( 2 )  pairs.  In  the  experimental  data  one 
can  observe  three  of  the  four  possible  crystal  orientations  of 
these  pairs  of  substitutionals. 

An  estimate  of  the  number  of  neighboring  pairs  of  Mn 
substitutionals  based  on  a  random  distribution  of  Mn  on  the 
Ga  sublattice  is  much  lower  than  the  number  of  compound 
features,  in  particular  those  labeled  “C,”  seen  in  Fig.  1.  One 
possible  explanation  is  positional  correlation  among  the  sub¬ 
stitutional  dopants  as  suggested  by  the  strong  short-range 
attraction  between  substitutional  Mn  found  in  a  recent  theo¬ 
retical  study.26  Moreover,  such  clusters  have  recently  been 
invoked  in  Monte  Carlo  calculations  to  investigate  the  phase 
diagram  of  Mn-doped  GaAs  with  quantitative  estimates  of 
the  ferromagnetic  ordering  temperature  in  agreement  with 
experiment.27'28  Nevertheless,  it  is  unlikely  that  a  description 
based  entirely  on  substitutionals  and  clusters  of  substitution¬ 
als  can  describe  the  frequency  and  shapes  of  all  the  features 
in  Fig.  1.  This  suggests  that  other  Mn  complexes,  in  particu¬ 
lar  those  involving  interstitial  Mn,  may  help  explain  the  re¬ 
maining  features. 


V.  DOPANT  COMPLEXES 
A.  Total-energy  considerations 

It  has  recently  been  shown  that  although  the  formation 
energy  of  interstitial  Mn  is  much  higher  than  that  of  substi¬ 
tutional  Mn  in  bulk  GaAs,  the  presence  of  the  (001)  growth 
surface  significantly  enhances  the  likelihood  of  Mn  occupy¬ 
ing  interstitial  sites.12  Thus  one  expects  that  both  substitu¬ 
tional  and  interstitial  Mn  will  be  present  in  typical  Mn-doped 
GaAs  samples,  as  shown  by  recent  measurements  of  the  Mn 
distribution  in  Mn-doped  GaAs.11  That  both  are  relevant  to 
the  observed  magnetic  and  transport  properties  of  Mn-doped 
GaAs  has  been  addressed  in  recent  theoretical  studies  of  the 
magnetic  interaction  in  nearest-neighbor  pairs  of  substitu¬ 
tional  and  interstitial  Mn,29  and  in  complexes  composed  of 
an  interstitial  Mn  trapped  between  two  neighboring 
substitutionals.’0  In  the  bulk,  the  barrier  for  diffusion  of  a 
charged  (  +  2e)  interstitial  Mn  is  —0.5  eV;  therefore  at  room 
temperature  we  expect  interstitial  Mn  to  rapidly  diffuse 
throughout  the  bulk  material  unless  it  becomes  kinetically 
trapped.  Since  substitutional  dopants  are  acceptors,  they  are 
a  natural  trap  for  interstitials,  electrostatically  binding  to 
them  in  a  number  of  different  physical  configurations.  Clus¬ 
tering  phenomena  for  compensating  dopants  such  as  these 
have  recently  been  predicted  using  Monte  Carlo 
techniques.31  On  the  other  hand,  since  interstitial  Mn  in  bulk 
is  a  double  donor  with  charge  of  +2e,  under  conditions  of 
filled-state  imaging  (for  which  the  tip  is  positively  biased 
relative  to  the  sample)  free  interstitial  Mn  will  be  repelled 
away  from  the  surface,  leaving  behind  only  those  interstitials 
that  are  bound  in  complexes.  To  determine  which  such 
interstitial-substitutional  complexes  are  most  likely  to  be 
present,  we  turn  to  total-energy  calculations. 
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FIG.  4.  Relative  total  energies  of  complexes  of  interstitial  and 
substitutional  Mn  near  the  GaAs  (110)  surface:  (a)  ,s(l  )iG‘‘(  I )  and 
s(l)tAs(l),  (b)  s(m)iGa(2)  and  s(m)iAs(  2)  (m=l,2),  (c) 

iGa(l)s(l)  and  iAs(l)s(l),  and  (d)  iG\2)s(l).  Simulated  STM 
images  are  provided  for  the  low-energy  configurations  in  the  shaded 
regions. 

To  simplify  the  discussion,  we  extend  the  notation  intro¬ 
duced  earlier  to  now  include  Mn  interstitials,  and  complexes 
of  interstitial  and  substitutional  Mn.  Isolated  interstitial  Mn 
in  layer  n,  coordinated  to  As  and  Ga,  will  be  referred  to  as 
iA\n)  and  iG\n),  respectively.  For  reference,  the  configu¬ 
ration  of  isolated  interstitial  dopants  in  the  first  two  layers  is 
shown  in  Fig.  2(b).  Complexes  composed  of  interstitial  and 
substitutional  Mn  will  be  referred  to  with  a  composite  nota¬ 
tion  that  takes  into  account  the  crystalline  orientation  of  the 
atoms  involved.  For  example,  since  the  GaAs  (110)  surface 
does  not  have  a  (001)  mirror-plane  symmetry,  there  are  two 
distinct  ways  of  forming  a  complex  of  interstitial  and  substi¬ 
tutional  Mn.  In  Fig.  2  the  surface  Ga  atoms  are  to  the  right  of 
the  surface  As  atoms  to  which  they  are  bonded,  so  that  we 
may  consider  complexes  where  the  interstitial  Mn  is  to  the 
left  or  right  of  the  neighboring  substitutional  dopant;  such 
complexes  will  be  denoted  as  iAs(n)s(m)  and  s(m)iAs(n), 
respectively.  Figures  2(c-e)  show  the  configuration  of  near¬ 
surface  complexes  we  have  considered  in  this  work.  Thus, 
for  example.  Fig.  2(c)  shows  the  complexes  s(2)iGd(2)  and 
s(2)/As(2),  which  we  discuss  below. 

Figure  4  shows  the  calculated  total  energies  of  Mn  com¬ 
plexes  relative  to  that  of  the  s(2)lA'(2)  dopant,  which  we 
identify  as  the  lowest-energy  configuration.  There  is  a  large 
variation  in  the  energies,  up  to  — 1.3  eV  above  the  energy  of 
s(2)/As(2).  We  also  note  that  for  complexes  involving  only 
Mn  atoms  on  the  surface  [Figs.  4(a)  and  4(c)]  the  lowest- 
energy  configuration  occurs  when  the  interstitial  site  is  near¬ 
est  the  substitutional  site.  Relative  to  these  nearest-neighbor 
complexes,  when  interstitial  Mn  is  displaced  by  lateral  trans¬ 
lations  in  the  (110)  direction,  the  energy  increases  dramati¬ 
cally;  the  next  lowest  energy  configurations  are  —0.3  eV 
higher  in  energy.  Thus,  diffusion  of  surface  interstitial  Mn 
from  high-energy  locations  (denoted  with  *  and  **)  to  low- 
energy  configurations  is  likely  to  form  nearest-neighbor 
complexes  on  the  surface. 
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FIG.  5.  Simulated  STM  images  for  Mn  substitutional-interstitial 
complexes.  The  locations  of  substitutional  and  interstitial  Mn  are 
denoted  with  +  and  X ,  respectively.  The  crystalline  orientation  is 
the  same  as  in  Fig.  2. 

It  is  difficult  to  estimate  on  energetic  grounds  alone,  the 
expected  number  of  such  complexes.  The  presence  of  both 
the  surface  and  the  STM  tip  create  additional  complications. 
For  example,  the  energy  of  interstitial  Mn  depends  strongly 
on  its  proximity  to  the  GaAs  (110)  surface;  an  interstitial  Mn 
in  the  interior  of  GaAs  far  from  the  surface  can  lower  its 
energy  as  much  as  ~  1  eV  by  moving  to  a  location  just  be¬ 
low  the  surface.32  This  suggests  that  interstitial  Mn  will  dif¬ 
fuse  from  the  bulklike  region  toward  the  surface;  during  the 
STM  scans,  these  “excess”  interstitials  will  be  observed  only 
if  they  are  bound  in  complexes.  Such  phenomena  are  diffi¬ 
cult  to  quantify  and  beyond  the  scope  of  this  work.  Instead 
we  choose  simply  to  report  simulated  STM  images  of  several 
low-energy  complexes  within  the  shaded  regions  of  Fig.  4. 

B.  Simulated  STM  Images 
1.  Substitutional-interstitial  complexes 

Figure  5  shows  the  simulated  images  for  the  low-energy 
configurations  of  substitutional-interstitial  complexes.  The 
simulated  image  of  the  lowest-energy  substitutional- 
interstitial  surface  complex  s(  1  )zGa(  1 )  is  shown  in  Fig.  5(a). 
There  is  a  reduction  in  intensity  of  the  As  surface  atom  to  the 
right  of  the  s(  1)  site  and  in  the  region  between  the  As  (001) 
planes  to  the  left  of  the  ,v ( 1 )  site.  The  simulated  images  of 
s(l)zGa(l)*  [Fig.  5(c)],  s(2)zAs(2)  [Fig.  5(d)],  and 
s(l)zAs(2)  [Fig.  5(f)]  are  all  similar,  but  the  reduction  in  the 
simulated  intensity  between  the  As  (001)  planes  to  the  left  of 
the  s(l)  site  is  less  apparent.  Among  these,  only  the 
s(2)zAs(2)  complex  has  a  (TlO)  mirror-plane  symmetry.  All 
four  of  these  simulated  images  are  generally  consistent  with 
the  characteristics  of  feature  E  in  the  experimental  data  of 
Fig.  1.  However,  since  the  s(2)zAs(2)  complex  is  much 
lower  in  energy  (~0.5  eV)  than  the  other  complexes  just 
described,  and  since  there  appears  to  be  only  two  such  fea- 


FIG.  6.  Simulated  STM  images  for  interstitial-substitutional 
complexes.  The  locations  of  substitutional  and  interstitial  Mn  are 
denoted  with  +  and  X ,  respectively.  The  crystalline  orientation  is 
the  same  as  in  Fig.  2. 

tures  in  Fig.  1,  we  suggest  that  feature  E  is  most  likely  a 
s(2)zAs(2)  complex.  If  this  is  correct,  then  the  observed  ap¬ 
parent  displacement  of  the  surface  As  atom  results  mainly 
from  an  electronic  effect  rather  than  a  real  atomic  displace¬ 
ment. 

The  simulated  image  of  the  s(l)/As(l)  complex  shown  in 
Fig.  5(b)  is  rather  different  than  those  for  the  complexes  just 
described.  This  simulated  image  obviously  has  a  ( 1 10) 
mirror-plane  symmetry  and  a  “butterfly”  shape  with  a  region 
of  lower  intensity  between  the  Mn  atoms.  Since  there  is  no 
(001)  mirror  plane,  the  left  and  right  wings  of  the  butterfly 
are  clearly  different,  with  the  contribution  from  the  ,v  ( 1 )  Mn 
more  intense  than  that  of  the  zAs(l).  Overall,  these  charac¬ 
teristics  are  similar  to  those  of  feature  C,  which  appears  to 
involve  four  As  surface  atoms  in  neighboring  (001)  planes. 
As  we  discuss  below,  there  are  other  energetically  competi¬ 
tive  complexes  with  similar  simulated  images,  so  that  we 
postpone  any  conclusions  regarding  feature  C  until  all  the 
relevant  complexes  have  been  discussed. 

The  simulated  image  of  the  s(2)zGa(2)  complex  shown  in 
Fig.  5(e)  does  not  have  a  (110)  mirror  plane  symmetry.  In 
this  configuration  there  is  an  oblate  region  of  intensity  along 
the  As  (001)  plane  and  a  reduction  in  the  surface  As  atom 
intensity  nearest  the  interstitial  site.  On  the  basis  of  the  avail¬ 
able  STM  data,  we  cannot  rule  out  the  possibility  that  some 
of  those  features  labeled  B  in  Fig.  1  result  from  s(2)zGa(2) 
complexes. 

2.  Interstitial-substitutional  complexes 

Figure  6  shows  the  simulated  images  for  the  low-energy 
configurations  of  interstitial-substitutional  complexes.  The 
simulated  images  of  the  zAs(l)s(l),  z'As(l)s(l)*,  and 
zCa(2)s(l)  complexes  are  all  very  similar  in  shape,  with  an 
increase  in  intensity  on  the  four  nearest  surface  As  atoms  and 
between  the  As  (001)  planes  to  the  left  of  the  s(l)  site  for 
the  first  two  and  to  the  right  of  the  s(l)  site  for  the  latter. 
Obviously,  the  zAs(l)i(l)  and  /Ga(2).v(  1 )  complexes  have  a 
(110)  mirror-plane  symmetry,  whereas  the  zAs(l)s(l)*  ap¬ 
pears  to  approximately  have  this  symmetry.  The  characteris¬ 
tics  of  these  images  are  similar  to  that  of  the  s(l)zAs(l) 
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complex  in  Fig.  5(b),  but  with  more  intensity  between  the 
relevant  As  (001)  planes  as  a  result  of  the  interstitial  location. 

The  simulated  image  of  the  /Ga(  1  )s(  1 )  complex  shown  in 
Fig.  6(c)  has  no  apparent  symmetry.  There  is  an  increase  in 
intensity  between  the  two  surface  As  atoms  nearest  the 
/Ga(l)  site,  whereas  the  s(  1)  Mn  is  only  very  weakly  visible. 
Features  such  as  this  are  not  apparent  in  the  experimental 
data,  although  XSTM  scans  of  lower  doped  samples  may 
help  clarify  whether  they  are  present. 

From  our  STM  data  and  the  theoretical  simulations,  we 
cannot  unambiguously  assign  feature  C  to  one  particular 
complex  among  the  plausible  alternatives:  s(l)tAs(l), 
tAs(l)s(l),  tAs(l)s(l)*,  ;Ga(2)s(  1 ),  i(l)s(l),  j(3)j(3), 
and  s(l)s(3).  Considering  the  large  number  of  such  features 
in  Fig.  1  and  the  comparable  energies  of  the  different  com¬ 
plexes,  it  is  likely  that  more  than  one  of  these  dopants  are 
present  in  the  scan  area  of  Fig.  1.  Measurements  on  samples 


with  much  smaller  Mn  concentration  are  ongoing  and  should 
be  able  to  address  these  issues  in  more  detail. 

VI.  SUMMARY 

In  summary,  we  have  combined  first-principles  calcula¬ 
tions  and  high-resolution  XSTM  measurements  to  character¬ 
ize  Mn-doped  GaAs.  XSTM  on  a  (110)  cleavage  plane  re¬ 
veals  five  reproducible  features  not  found  in  bulk  GaAs. 
Total-energy  calculations  were  used  to  screen  the  most  likely 
configuration  of  Mn  dopants  near  the  GaAs  (110)  surface. 
The  low-energy  configurations  were  then  further  scrutinized 
by  comparison  of  simulated  and  measured  STM  images. 
These  comparisons  reveal  that  there  are  predominantly  two 
types  of  Mn-related  dopants  in  the  sample:  (1)  those  involv¬ 
ing  only  Mn  substitutional(s),  including  both  isolated  and 
pairs  of  substitutionals  and  (2)  complexes  composed  of  sub¬ 
stitutional  and  interstitial  Mn. 
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